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Abstract 

It is widely recognized that global warming promotes soil organic carbon (SOC) decomposition, and soils thus emit more 
CO2 into the atmosphere because of the warming; however, the response of SOC decomposition to this warming in 
different soil textures Is unclear. This lack of knowledge limits our projection of SOC turnover and CO2 emission from soils 
after future warming. To investigate the CO2 emission from soils with different textures, we conducted a 1 07-day incubation 
experiment. The soils were sampled from temperate forest and grassland in northern China. The Incubation was conducted 
over three short-term cycles of changing temperature from 5°C to 30°C, with an interval of 5"C. Our results indicated that 
CO2 emissions from sand (>50 ^m), silt (2-50 ^im), and clay (<2 |xm) particles Increased exponentially with increasing 
temperature. The sand fractions emitted more CO2 (CO2-C per unit fractlon-C) than the slit and clay fractions In both forest 
and grassland soils. The temperature sensitivity of the CO2 emission from soil particles, which is expressed as Qio, decreased 
In the order clay>sllt>sand. Our study also found that nitrogen availability in the soil facilitated the temperature 
dependence of SOC decomposition. A further analysis of the incubation data indicated a power-law decrease of Qio with 
Increasing temperature. Our results suggested that the decomposition of organic carbon In fine-textured soils that are rich 
In clay or silt could be more sensitive to warming than those in coarse sandy soils and that SOC might be more vulnerable in 
boreal and temperate regions than In subtropical and tropical regions under future warming. 

Citation: Ding F, Huang Y, Sun W, Jiang G, Chen Y (2014) Decomposition of Organic Carbon in Fine Soil Particles Is Likely IVlore Sensitive to Warming than in 
Coarse Particles: An Incubation Study with Temperate Grassland and Forest Soils in Northern China. PLoS ONE 9(4): e95348. doi:10.1371/journal.pone.0095348 

Editor: Dafeng Hui, Tennessee State University, United States of America 

Received December 3, 2013; Accepted IVlarch 26, 2014; Published April 15, 2014 

Copyright: © 2014 Ding et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was jointly supported by the National Natural Science Foundation of China (Grant No. 41075108), the Ministry of Science and Technology of 
China (2010CB950604) and the CAS Strategic Priority Research Program (Grant No. XDA05050507). The funders had no role in study design, data collection and 
analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 
* E-mail: sunwj@ibcas.ac.cn 



Introduction 

As the largest terrestrial carbon pool, soil plays a vital role in the 
global carbon cycle [1]. A small change in the magnitude of soil 
carbon stocks could result in a large impact on the atmospheric 
enrichment of CO2 [2]. Soil organic carbon (SOC) loss via soil 
heterotrophic respiration (i.e., SOC decomposition) is one of the 
largest components of global CO2 flux from the terrestrial 
ecosystem to the atmosphere [3,4]. It is widely recognized that 
temperature has a profound influence on SOC decomposition [5] 
and that soils are expected to release CO2 at faster rates with rising 
temperature [3]. Therefore, understanding the temperature 
sensitivity of SOC decomposition is critical for predicting changes 
in global soil respiration, changes in soil carbon stocks to future 
warming, and especially, the feedback of soils on climate change 
[2,6-8]. 

The temperature sensitivity of SOC decomposition, which is 
commonly expressed as Qjo, refers to the change in decomposition 
rate with a 10°C increase in temperature. A recent study 
constructed a database of worldwide soil respiration observations 
matched with high-resolution historical climate data and estimated 



that global soil respiration had a Qjo of 1.5 [3]. Nevertheless, soil 
respiration experiments have unanimously indicated that Qjo is 
spatially heterogeneous on a global scale [9,10]. Thus, to 
accurately predict changes in global soil carbon stocks to future 
warming, knowing which factors determine the temperature 
dependence of SOC decomposition on large scales is critical [6]. 
It is known that the Qjo of soil heterotrophic respiration varies 
with temperature [5,11], soil moisture [12,13], and the quantity 
and quality of SOC [14,15]. However, it is not clear whether SOC 
decomposition has different warming responses in different soil 
textures. It is most likely so because the SOC stock in fine-textured 
soils is less vulnerable to disturbances, e.g., tillage, than in coarse- 
textured soils [16,17]. 

Such different vulnerabilities are assumed to result from fine- 
textured soils having greater fractions of silt-sized or clay-sized 
particles, whereas coarse-textured soils have more sand-sized 
particles; these primary particles contain different mineral 
compositions and provide various afTmities to organic matter 
[18]. Clay and silt particles (e.g., sesquioxides and layer silicates) 
provide large specific surface areas and numerous reactive sites at 
which SOC can be sorbed by strong ligand exchange and 
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polyvalent cation bridges, whereas the sand particles, which is 
dominated by quartz particles, exhibits only weak bonding 
affinities to SOC [19,20]. The bonding of SOC to minerals is 
recognized to prevent decomposition [21,22] and may therefore 
decrease the temperature sensitivity of decomposition [23]. Thus, 
the decomposition of organic carbon associated with sand, sUt, and 
clay fractions likely respond to warming differently. The temper- 
ature sensitivity of SOC decomposition from different particle size 
fractions is intimately linked to the response of CO2 emission from 
different soil textures to future warming and the emission's 
feedback on climate change. 

Different soil particle size fractions possess different sizes of 
carbon pools [20]. Clay and silt have large specific surface areas 
and can sorb abundant SOC [24]. For example, in temperate 
arable soils, 50-75% of the total SOC is associated with clay, 
whereas silt accounts for 20-40% and sand accounts for less than 
10% [18]. Additionally, the SOC associated with the clay and sik 
fractions is assumed to be stable relative to that associated with the 
sand fraction [21]. In this context, the long-term stabihzation of 
soil carbon pools under climate warming mainly depends on the 
warming response of the carbon that is associated with the clay 
and sUt fractions. 

However, to our knowledge, there were few studies who 
explored difference in the warming responses of organic carbon 
associated with these primary particles. In this study, we 
conducted a 107-day incubation experiment to investigate the 
temperature sensitivity of SOC decomposition in sand (>50 |J,m), 
silt (2-50 nm), and clay (<2 |J.m) particles from grassland and 
forest soils. This study aims to answer the following questions: 1) Is 
there a difference in the Qjo of SOC decomposition among 
different size soil particles and between the two different 
ecosystems? 2) Which factors determine the Qjo? A further 
objective is to better understand the features of CO2 emission and 
their future trends with global warming for different textured soils. 

Materials and Methods 

Sites and sampling 

Soil samples were collected from two sites at the typical steppe 
grassland and temperate forest ecosystems in northern China. The 
grassland site was at the Duolun Restoration Ecology Experimen- 
tation and Demonstration Station (42°02'N, 116°17'E, 1324 m 
a.s.l.) and the forest site was at the Beijing Forest Ecosystem 
Research Station (39°58'N, 115°26'E, 1252 m a.s.l.) at the 
Institute of Botany of the Chinese Academy of Sciences. The 
grassland is permanent grassland situated in the temperate 
continental semi-arid climate zone. The dominant species are 
Stipa baicalensis, Thalictrum petabideum, and Artemisia Jrigida. The 
mean annual temperature (MAT) is 1.6°C, and the mean annual 
precipitation (MAP) is 386 mm, with 100 frost-free days. The soil 
is sandy Haplic Calcisols (FAO) composed of 61.2% sand, 23.0% 
silt, and 15.8'/'o clay, and the pH is 7.0. The forest is an artifitdal 
forest (with age >40 years) with a dominant pinus tabuliformis 
Carriere specie and the climate is warm temperate humid monsoon 
climate zone. The MAT and MAP are 5°C and 650 mm, 
respectively, with 160 frost-free days. The soil is a clay loam that 
belongs to Eutric Cambisol (FAO) and is composed of 41.0% sand, 
42.5% silt, and 16.5% clay, and the pH is 7.0. 

At each site, 40 cores of mineral soils at 0-20 cm depths were 
randomly collected using an auger with 4 cm in diameter and 
were composited into one soil sample. A small fraction of the soils 
(approximately 1 kg) was immediately placed in a sterile plastic 
bag on ice when transported to the laboratory and stored at 4°C 
for inoculum preparation. The remaining soils were placed into 



cloth bags and air-dried at room temperature. The dried soils were 
finely crumbed and sieved through a 2-mm mesh to remove visibly 
identifiable plant debris and gravel. Fine roots were removed 
through static electricity sorption. 

Particle-size fractionation 

Twenty grams of dried soil were dispersed in 100 ml of 
deionized water by the ultrasonic probe of a JY92-II Ultrasonic 
Homogenizer (6 mm diameter probe, 300 W'). The ultrasonic time 
was chosen as 3 min and 8.5 min for grassland and forest soil, 
respectively, after calibrating against conventional soil texture 
analyses [25] . The dispersed suspensions were separated into three 
size fractions: sand (>50 \ixn), silt (2-50 ^im), and clay (<2 nm). 
The sand fraction was separated using distilled water to wash all of 
the fine particles in the soil suspension through a 50-Hm sieve until 
the elutriate solution that exited the sieve was clear. The elutriate 
was divided into silt and clay fractions by repeated centrifugation 
using an Avanti J-25 Beckman Centrifuge [26]. The clay fraction 
was concentrated from the large volumes of suspensions by 
flocculation with CaClg and centrifugation [27]. The three size 
fractions were oven-dried at 40°C [28], gently crushed, weighed, 
and then stored for incubation and future analysis. 

Laboratory incubation 

For each soil particle, a 10-g dried sample was placed into a 
150-ml brown flask. Then, 2 ml of inoculum were added to each 
sample, and the inoculum solution was prepared by shaking 1 00 g 
of fresh soil with 1 000 ml of deionized water and leaving to stand 
overnight [28]. To improve soil aeration and make the inoculum 
disperse uniformly, the sample plus inoculum was loosened, 
stirred, and leveled slightly. The soil moisture was adjusted to 60% 
water-filled pore space (WFPS) [29]. The incubation lasted 107 
days (Figure 1). The soil samples were incubated at room 
temperature (ca. 20°C) when there was no measurement. The 
samples were regularly corrected for water loss by adding 
deionized water. Four replicates were conducted in the experi- 
ment. Empty flasks were set as blank. The volume of head space 
(flask volume minus soil solid) was accurately measured after 
complementing the experiment by the method of water replace- 
ment [30]. 

Following Fang et al. [31] and Xu et al. [32], soil respiration 
rates under changing temperature were measured in sequence. 
Three temperature-change cycles were started at the 9*, 47'*', and 
lOO**" days after the beginning of the incubation, and each cycle 
lasted 8 days (Figure 1). The first temperature-change cycle began 
after 8 days of incubation to stabilize the microbial activity to 
avoid an undesired microbial peak [33]. Within each cycle, the 
temperature was continuously increased from 5°C to 30°C in 
intervals of 5°C and then gradually decreased from 30°C to 5°C 
with the same interval. After each temperature was changed to the 
next, the incubation flasks were kept open for 3 hours to allow the 
soil samples to adapt to the changed temp(;rature [30]; the flasks 
were then sealed using butyl rubber septa and flushed with 
ambient air for one minute using a pump with an "inlet-outlet" 
needle inserted through the septa to produce the same initial gas 
conditions among all flasks. After specific times of enclosure (20, 
15, 9, 6, 4, and 3 hours for 5, 10, 15, 20, 25, and 30°C, 
respectively), 5 ml of gas was taken through a gas-tight injection 
syringe from the headspace of the flask. The CO2 concentrations 
were analyzed with a gas chromatograph (Agilent 7890A, USA). 

Determination of soil physicochemical properties 

The soil texture was measured using the hydrometer method 
[34] . The SOC and soil total nitrogen (STN) concentrations were 
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Figure 1. Schematic diagram of tKie time sequence for tlie tliree temperature-cliange cycles. RT signifies that the soil samples were 
incubated at room temperature. 
doi:l 0.1 371 /journal.pone.0095348.g001 



determined using the K2Cr207-H2S04 wet oxidation method [35] 
and the Kjeldahl method [36], respectively. The microbial 

biomass carbon (MBC) and microbial biomass nitrogen (MBN) 
in the .soil samples at the end of incubation were measured using 
chloroform fumigation-extraction (K2SO4) [37] with a Multi N/C 
3100-TOC/TN Analyzer (Analytikjena, Germany). The K2SO4- 
extractable organic carbon and nitrogen from the non-fumigated 
soil samples wcr(^ considered as the dissolved organic carbon 
(DOC) and total dissolved nitrogen (TDN) [38,39]. 

Calculations and statistical analysis 

Similar to Lang et al. [40], the CO2 emission rate was calculated 
using the following equation: 

273 12 
^=^'' (273 + 71 ^^'^-"'^^^^^^44/^^^^^/^ 

where R is the CO2 emission rate ()j.g CO2-C g~' OC li '); Fis 
the flask headspace volume (L); T is the incubation temperature 
(°C); C, and Q, are the CO2 concentrations in the sample and 
blank flasks (ppm), respectively; p is the standard state CO2 density 
(g L~'); li'is the dry weight of the soil particles (kg); c is the SOC 
concentration of the soil particles (g kg '); and / is the time 
enclosed (h). The final CO2 emission rate at a given temperature 
was calculated as the average of the values between the increasing 
temperature and decreasing temperature periods [31,32]. To 
illustrate the temperature sensitivity of CO2 emission from soil 
particles, we used two methods to calculate Qjo. One method 
involves the commonly used first-order exponential model: 



R = Ae^ 



(2) 



where R is the CO2 emission rate (|ig CO2-C g~' SOC h~'), Tis 
the temperature (°C), and A and B are model parameters that are 
to be determined. The parameter A is a carbon quality index, 
which refers to the availability and lability of the carbon substrates 
[15,41]. The Qjo value was then computed as following: 



Qio = e" 



(3) 



The Qjo calculated by Eqn. 3 was a constant value in 
temperature range of 5-30°C. To examine the Qjo variation 
with temperature, the Qjo values were also calculated separately 
for each temperature interval of 5 degrees (e.g., 5-1 0°C) and 10 
degrees (e.g., 5-1 5°C) based on the following equation [11]: 



(4) 



where R2 and Ri are the CO2 emission rates (ng CO2-C g SOC 
h ') measured at temperatures T2 (high temperature) and Ti (low 
temperature), respectively. 



All statistical analysis was conducted using SPSS 13.0 (SPSS 
Inc., 2004). The means were compared using the Duncan's 
multiple range test (DMRT). 

Results 

Physicochemical characteristics of soils 

The SOC and STN concentrations of the bulk soil were very 
similar between the grassland and forest (Table 1). The organic 
carbon concentrations in grassland and forest soil particles both 
decreased with increasing particle size, from 55 g kg~' and 73 g 
kg~ in the clay particle (<2 |im) to 9 and 13 g kg~ in the sand 
particles (>50 |xm), respectively. Similar trends were observed for 
the total nitrogen concentrations. The C:N ratios were similar 
between the clay and silt particles, which were lower than those of 
the bulk soil and the sand particles. The total percentages of the 
carbon massc's associated with clay and silt particles in the total 
SOC were more than 60% for both grassland and forest soils, and 
the percentages of their nitrogen masses in the total STN were 
both approximately 75%. 

Relationship between soil organic carbon decomposition 
rate and temperature 

The SOC decomposition rate from soil particle size fractions, 
which is expressed as the CO2-C emission rate per unit of fraction- 
C, had an exponential relationship with the incubation temper- 
ature in all measurements (Figure 2). The sand fractions, which 
had a lower organic carbon content (Table 1), exhibited greater 
SOC decomposition rates than the silt and clay fractions, which 
latter two were similar to one another. The decomposition rates in 
all the soil fractions decreased continuously as incubation time 
progressed (from an earlier temperature-change cycle to a later 
one), except for the sand fraction from the second to the third 
cycle. As a whole, the decomposition rates from grassland soil 
particles were greater than those from the corresponding forest soil 
particles under same condition (Figure 2). 

The parameter A in the fitted equation (Eqn. 2) for the 
relationship between the decomposition rate and temperature is an 
index for carbon availability and lability. The value of parameter 
A was greater for the sand fraction than for the silt and clay 
fractions in both the grassland and forest soils (Figure 2). 
Furthermore, the values of parameter A were approximately 
erjual between the two land uses for the clay and silt fractions 
during each temperature change cycle but were greater in the 
grassland than in the forest for the sand fraction (Figure 2). 

Temperature sensitivity of soil organic carbon 
decomposition for different soil particle-size fractions 

The Qio values calculated using the exponential model (Eqn. 3) 
were significandy different among the soil particle size fractions 
during each temperature-change cycle; the clay fraction exhibited 
the greatest value, followed by the silt and sand fractions (Table 2). 
The calculated mean Qio values of the three cycles were 
2.68±0.16 and 2.43±0.16 for the clay fraction, 2.45±0.06 and 
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Table 1. Physicochemical composition of the soils. 





Land uses and particle size fractions 


Fraction ' 


(%) Organic carbon 




Total nitrogen 




C:N ratio 






Concentration (g 
kg ') 


Share of the total Concentration (g 
(%) kg ^ 


Share of the total 

(%) 




Grassland 


Clay (<2 fim} 


11.2 (1.6) 


55.3 


28.5 


5.8 


31.9 


9.5 


Silt (2-50 |im) 


28.0 (1.3) 


30.9 


40.0 


3.1 


43.2 


9.9 


Sand (>50 urn) 


57.4 (1.0) 


8.9 


23.6 


0.7 


18.7 


13.4 


Bulk soil 




21.7 




2.0 




10.7 


Recovery (%) 


96.6 (1.1) 




92.1 




93.9 




Forest 


Clay (<2 |im) 


18.1 (2.8) 


72.6 


44.6 


6.1 


51.2 


12.0 


Silt (2-50 |jm) 


35.1 (2.7) 


17.4 


20.8 


1.5 


24.2 


11.8 


Sand (>50 ^lm) 


43.9 (1.2) 


13.0 


19.4 


0.7 


14.1 


18.9 


Bulk soil 




29.5 




2.1 




13.8 


Recovery (%) 


97.1 (1.0) 




84.8 




89.5 





^Values in the parentheses represent SD. 
doi:l 0.1 371 /journal.pone.0095348.t001 



2.22±0.10 for the sUt, and 2.33±0.08 and 2.09±0.09 for the sand, 
in the grassland and forest soils, respectively. The mean Qjo values 
of the grassland soil particles were significandy greater than those 
that of the corresponding forest soil particles (Table 2). As the 
incubation time progressed (from an earlier temperature-change 
cycle to a later one), Qjo increased for the grassland clay and silt 
fractions (Table 2). 

The Qjo values decreased in power-law with increasing 
temperature for all soil particles when calculating using Eqn. 4 
during each individual temperature interval (Figure 3). The 
average Qjo values of all soil particle fractions were 4.22±0.14 
during the minimum temperature interval and 1.81±0.17 during 



the maximum temperature interval (Figure 3). The grassland soil 
particles had greater Qjo than the corresponding forest soil 
particles during the temperature interval of 5-15°C, whereas the 
Qjo values for the two land uses tended to approach one other 
during 15-30°C (Figure 3). Nevertheless, the grassland soil 
particles exhibited a more rapid decline in Qjo with increasing 
temperature (i.e., greater absolute values of the power-law 
exponents in Figure 3) than did the forest soil particles. 
Furthermore, the Qjo values calculated using Eqn. 4 followed as 
clay>silt>sand (Figure 3), which was consistent with the pattern 
in Qjo calculated using the exponential model (Eqn. 3) (Table 2). 



*/ -o 7=0.087x 

■ y3ii,=2.0e'''"88x 
■y,„,=2.0e°°^'"< 



R==0.98 y^O.001 
R==0,98 P<0.00^ 
R==0.98 /OcO.OOl 
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ysii,= 1,0e"™"'' R==0.97 P<0.001 
y^,i3y=0.6e°°i''»^ R^=0.98 /O<0.001 
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Figure 2. Relationship between C02 emission rates from particle-size fractions and incubation temperature for grassland (A, B, C) 
and forest soils (D, E, F). The panels from left to right represent temperature-change cycles in the order 1, 2, 3. The vertical bars denote the 
standard deviations (n=4). 
doi:10.1371/journal.pone.0095348.g002 
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Table 2. The Qio" values {±SD) of SOC decomposition for different particle-size fractions in grassland and forest soils {n = 4). 



Land uses 


Temperature-change cycles 


Particle size fractions 










Clay «2 (im) 


Silt (2-50 iim) 


Sand (>50 )im) 


Grassland 


1 


2.56±0.07a'' 


2.41±0.a5b 


2.39±0.a4b 




2 


2.60±0.04a 


2.45±0.06b 


2.34±0.09c 




3 


2. 88 ± 0.04a 


2.51±0.02b 


2.27±0.07c 




Mean* 


2.68±0.16aA* 


2.45±0.06bA 


2.33±0.08cA 


Forest 


1 


2.45±0.07a 


2.32±0.08b 


2.19±0.04c 




2 


2.24±0.01a 


2.13±0.08b 


2.08±0.03b 




3 


2.59±0.04a 


2.23±0.a5b 


2.00±0.05c 




Mean* 


2.43±0.16aB 


2.22±0.10bB 


2.09±0.09cB 



*Qio value was calculated using Eqn. 3 with the parameters for the fitted curves in Figure 2. 
*The mean was calculated including each replicate among the three temperature-change cycles (n = 12). 
^'Different lowercase letters in the same row indicate significance differences (p<0.05} among soil fractions. 
^Different capital letters in the same column indicate significance differences (p<0.05) between two land uses. 
doi:l 0.1 371 /journal.pone.0095348.t002 



Dependence of Qio on soil biochemical properties 

Generally, there were large differences among various particle- 
size fractions in most soU biological and chemical properties at the 
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Figure 3. Dependence of Q10 on incubation temperature for 
(A) clay, (B) silt, and (C) sand. The QIO for each temperature interval 
of 5 degrees (e.g., 5-10°C) and 10 degrees (e.g., 5-15'C) was calculated 
using Eqn. 4. The average temperature for each interval was taken as 
representative temperature. Data from all the three temperature- 
change cycles were included. Vertical bars denote the standard errors 
(n = 12). 

doi:10.1371/journal.pone.0095348.g003 



end of incubation (Table 3). Specifically, MBC, MBN, DOC, and 
TDN were greatest in the clay fraction and were generally of the 
same order of magnitude as Qjo (Table 2). However, the ratios of 
MBC/OC and MBN/TN shifted to be greatest in the sand 
fraction. The microbial C/N ratios were not significantly different 
among the soil particles, but they were generally greater than the 
soil particles' C/N ratios (except for the forest sand particles; see 
Table 1). Additionally, the grassland soil particles had greater 
MBC/ OC and less qC02 values than the corresponding forest soil 
particles. 

The soil biological and chemical properties at the end of 
incubation represented the soil conditions during the 3"* 
temperature-change cycle; thus we analyzed the relationship 
between the Qjq in the 3'^'' cycle and these soil biological and 
chemical properties. Table 4 presents the Pearson correlation 
coefficients of Qjo during an intermediate temperature range of 
15-25°C with various soil biochemical properties. Clearly, the Qjq 
values were highly correlated with MBC and TDN. The Qjo had 
significant linear relationships with MBC (Figure 4A) and TDN 
(Figure 4B) for the 5-1 5°C and 20-30°C temperature intervals 
and also for other temperature intervals (not shown). It should be 
noted that the slopes of regression line were different between the 
temperature ranges of 5-15°C and 20-30°C. Over an increasing 
temperature interval of 10°C, the slopes with low baseline levels of 
5°C were steeper than those with high baseline levels of 20°C. 
Note that the datasets were weU-mrxed for the grassland and 
forest, thereby suggesting that these two land use soils had similar 
features in terms of the dependence of Qjo on MBC and TDN. 

Discussion 

Carbon stability and microbial activity among different 
size fractions 

The carbon associated with different particle-size fractions has 
different availability and lability [42]. The sand fraction had 
greater organic carbon availability and lability (greater parameter 
A), and thus lower carbon stability, than did the silt and clay 
fractions (Fig. 2). This pattern of SOC stability among the different 
size fractions was consistent with previous studies [43,44]. The 
difference in SOC stability among soil primary particles is the 
reason why coarse sandy soils that are rich in sand particles are less 
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stable and exhibit faster SOC turnover rates than fine-textured 
soils that are rich in clay and silt particles [45] . 

Land-use type can affect the SOC stability by inputting carbon 
into the soil via litter or roots with different chemical recalcitrance 
[18,46]. The sand fraction had greater carbon stability (smaller 
parameter A) in the forest than in the grassland soil (Fig. 2). A 
possible reason is that the organic matter in the sand fraction is 
primarily composed of partially decomposed plant litter [18]; 
forest litter has higher lignin concentration and higher C:N ratio 
than grassland litter and is thus usually more recalcitrant [47] . In 
contrast, the carbon stability (parameter A) in clay or sUt fractions 
was similar between grassland and forest during each temperature- 
change cycle (Fig. 2), thereby suggesting that the organic carbon 
stabilities in clay and silt fractions were not affected by land-cover 
type. The carbon stability in these fine fractions is primarily 
related to organo-mineral association (i.e., chemical protection) 
[48]. Overall, the SOC stability in the sand fraction is primarily 
determined by the chemical recalcitrance of organic matter but is 
more strongly determined by the association of SOC with minerals 
in the clay or silt fractions. This pattern was consistent with 
different determinants of organic matter-derived nitrogen preser- 
vation between the sandy soils and clay soils [49] . 

Microbial activities and distribution may differ among different 
size fractions because these fractions provide distinct microhabitats 
and various organic carbon qualities and quantities [50-52] . Clay 
fraction had the greatest organic carbon concentration (Table 1) 
and had the maximum amount of MBC (Table 3), thereby 
indicating that the greatest biological activity was in clay fractions 
[18]. However, the clay fraction had the least MBC/OC ratio 
(Table 3), an eco-physiological parameter that refers to the 
microbial efficiency in using carbon substrate [53]. The microbes 
in the clay fraction had the lowest carbon use efficiency, which is 
likely related to low carbon availability due to strong organo- 
mineral association. Moreover, the grassland soil particles had 
greater MBC/OC and lower qC02, i.e., microbial community 
respiration per unit microbial biomass [53], than the correspond- 
ing forest particles (Table 3), thereby suggesting that the grassland 
soil contains a microbial community that is more efficient at using 
carbon energy. 

Effect of climatic and soil parameters on Qio 

First, temperature can determine the Qjo of SOC decomposi- 
tion. The lower Qjo of soil respiration at higher temperatures was 
not only assumed by the Arrhenius equation [6] but also 
confirmed by two comprehensive analyses of incubation studies 
of soils from all types of terrestrial ecosystems [11,54]. Similarly, 
we found that the Qjo values were negatively correlated with the 
incubation temperature (Figure 3). The relationship was fitted by a 
power function: y = ax^. Clearly, higher values of 0( produce 
higher values of Qjq. Otherwise, with larger absolute values of P, 
Qjo decreases more rapidly with increasing temperature. The 
parameters a and (5 in the Qjo-T curves were found to be 
significantiy correlated with the MBC and soil C:N ratio (Figure 5). 
An increase in the MBC concentration would lead to ot and Qjo 
rising (Figure 5A), which is consistent with Figure 4A. Significant 
inverse relationships were found between the parameter a and the 
soil C:N ratio (Figure 5B), thereby suggesting that higher soil C:N 
ratios was linked with higher tolerance to increasing temperature. 
A more rapid rate of Qjo decrease with increasing temperature 
(greater absolute value of fi) in grassland versus forest soil particles 
could be related to lower C:N ratios (Table 1) or greater MBC in 
grassland soil particles (Table 3). It should be mentioned that no 
significant relationship was found between the values of Qt and P 
and the soil properties of MBN, DOC, and DTN. 
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Table 4. Pearson correlation coefficient for Q,o {15-25°C) in the third temperature-change cycle and soil biological and chemical 
properties. 







Q,o (15-25°C) • 


MBC 


MBN 


DOC 


TDN 


Micro C/N 


DOC/TDN 


MBC/OC 


MBN/TN 


qCO^ 


Qio(15-25°C) 1 


MBC 


0.87"" 


1 


















MBN 


0.52" 


0.44* 


1 
















DOC 


0.58"" 


0.32 


0.38 


1 














TDN 


0.90""" 


0.80*** 


0.45" 


0.81""" 


1 












Micro C/N 


0.51" 


0.59** 


-0.29 


0.25 


0.54"" 


1 










DOC/TDN 


-0.85*"" 


-0.73*** 


-0.42 


-0.20 


-0.59"* 


-0.45" 


1 








MBC/OC 


-0.02 


0.19 


-0.09 


0.65"" 


0.29 


0.13 


0.19 


1 






MBN/TN 


-0.43* 


-0.29 


0.04 


-0.59"" 


-0.54"* 


-0.44" 


0.19 


0.71"* 


1 




qC02 


-0.66""" 


-0.70*"" 


-0.24 


0.18 


-0.35 


-0.37 


0.85""" 


-0.54* 


-0.05 


1 



^Qio (15 — 25'C) was calculated using Eqn. 4 for each replicate. 

*, **, ***Significant at probability levels of 0.05, 0.01, and 0.001, respectively. 

doi:l 0.1 371 /journal.pone.0095348.t004 



The negative relationship between Qjq of SOC decomposition 
and temperature has two implications. First, in a given region, 
SOC decomposition should be more sensitive to warming during 
winter than warming during summer. This trend was observed in 
the Alaskan tundra, with an increase of greater than 50% in CO2 
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Figure 4. Dependence of Q10 in the third temperature-change 
cycle on microbial biomass carbon (A) and total dissolved 
nitrogen (B). The Q10 (5-15°C) and Q10 (20-30 'Q were calculated 
using Eqn. 4 for each replicate. 
doi:10.1371/journal.pone.0095348.g004 



fluxes during the snow-covered period by winter warming and a 
20% increase during the growing season by summer warming 
[55]. Nevertheless, considering that the initial soil respiration rate 
is much less in winter than in summer, the absolute magnitude of 
the CO2 emission increase from warming could be still less in 
winter than in summer. Second, the SOC stock in boreal and 
temperate regions (low MAT), especially peat lands and perma- 
frost, should be more vulnerable to warming than that in 
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Figure 5. Dependence of parameters in the Q10— T curve on 
the soil microbial biomass carbon (A) and soil C: N ratio (B). 

Vertical bars denote the 95% CI of the estimated parameters in Fig. 3. 
For each parameter oc and p, six points represent the three particles of 
clay, silt, and sand from the grassland and forest soils. Horizontal bars 
represent the standard errors (n = 4). 
doi:1 0.1 371/journal.pone.0095348.g005 
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subtropical and tropical regions (high MAT) [6]. Moreover, peat 
lands and permafrost contain more than 50% of global soil carbon 
stock [56] and have been experiencing remarkable warming [57]. 
Consequendy, peat lands and permafrost are expected to 
experience largest potential soil carbon loss by 2 1 00 due to global 
warming [6]. 

Second, soil texture could have an impact on Qjo of SOC 
decomposition. The Qjo decreased in the order clay>silt>sand 
seemingly regardless of land-cover types, soil types, and climatic 
factors (Table 2). Consequently, SOC decomposition in fine- 
textured soils (rich in clay and silt content) could be more sensitive 
to global warming than that in coarse-textured soils (rich in sand 
content). Moreover, fme-textured soils generally have greater 
capacities to preserve SOC [58] and thus contain more SOC than 
do coarse-textured soils [1,59]. Therefore, current global model 
predictions that use constant Qjo values for different soil textures 
may predict an underestimated response of the soil carbon cycle to 
climate warming. 

Third, the soil nitrogen availability, indicated by TDN, could 
also affect the Qjo values, and greater Qjo values were related to 
higher TDN content regardless of land use type (Figure 4B). This 
result agrees with other incubation studies: Chang et al. [60] 
demonstrated the existence of a positive relationship between Qjo 
and the soil total nitrogen content, and Cusack et al. [61] observed 
that nitrogen fertilization significantly increased the Qjo of slow 
SOC pools. This relationship between Qjo and nitrogen 
availability could result from the nitrogen availability driving the 
activities of the soil enzymes that acquire carbon [62]. Low 
nitrogen availability may inhibit enzyme activity and thus 
constrain the warming response of SOC decomposition. This 
close relation of Qjo with soil nitrogen availability suggests that 
atmospheric nitrogen deposition could stimulate the positive 
feedback of soil respiration to climate warming. 

Qio for labile carbon pools vs. Qio for recalcitrant carbon 
pools 

The question of whether the decomposition of recalcitrant and 
labile carbon has different temperature dependencies is important 
because most experimental work has been done on the temper- 
ature scnsiti\'ity of the decomposition of more labile fractions, yet 
the long-term effect of warming on soil carbon storage should be 
primarily determined by the temperature dependence of the more 
recalcitrant carbon pool [7]. Kinetic theory and the Arrhenius 
equation both suppose that recalcitrant carbon pools are more 
warming-sensitive than labile carbon pools [6,6,3]. Many labora- 
tory incubation studies have supported this supposition [32,64,65]. 
Our study further affirmed this result by showing that the 
decomposition of the carbon associated with the clay and silt 
fractions (recalcitrant carbon pool), with slower turnover rates 
(Figure 2), yet had higher Qjo than that with the sand fraction 
(labile carbon pool) (Table 2). Furthermore, the organic carbon 
associated with the clay and silt fractions usually accounts for a 
large portion of the total SOC, i.e., more than 6Q"/o in this study 
(Table 1). Therefore, long-term response of th(- SOC stock to 
cUmate warming would be more severe than the current 
predictions that assume the same Qjo for labile and recalcitrant 
carbon pools. 

However, the supposition that recalcitrant carbon pools have 
greater warming dependence has not obtained wide acceptance 
[31,66]. Based on this supposition, Qjo in the incubation studies is 
expected to rise with prolonged incubation time because SOC is 



generally composed of a large proportion of labile carbon during 
the early incubation period, yet of a small proportion of labile 
carbon during the late incubation period due to depletion of labile 
carbon [67]. And yet, Fang et al. [31] and Reichstein et al. [66] 
found that Qjo was stable and independent of the incubation time. 
Similarly, our study only observed increasing trends in Qjo for the 
grassland clay and silt fractions from the T' temperature-change 
cycle to the 2"'' and the 3'^'' (Table 2) and did not found such trends 
in Qjo for the forest soil fractions or the grassland sand fractions. 
The inconsistency between these observations and the expectation 
derived from the above supposition could be explained by that Qjo 
having a significant positive relation with MBC and lower MBC 
values resulting in a smaller Qjq (Figure 4A); MBC generally 
declines with prolonged incubation because of the depletion of 
labile carbon [31,68]. Consequently, as incubation is prolonged, 
increasing trend of Qjo due to increased recalcitrant carbon 
proportion could be counteracted by decreasing trend of Qjo as a 
result of declined MBC [67]. If the MBC did not decrease sharply 
from early to late incubation stage, the Qjo would have increased 
as incubation progressed in Fang et al. [31] and Reichstein et al. 
[66], as well as in our study. However, the mechanism for how Qjo 
is affected by MBC is not clear. Assuming that the microbial 
community composition does not change and that every single 
microbe maintains the same warming response with prolonged 
incubation time, a decline in MBC alone should not induce lower 
Qjo. Possible mechanisms could be either microbial community 
shift to species with low thermal response [69] or microbe 
acclimatization to temperature increase [70]. On the other hand, 
soil carbon recalcitrance did not certainly increase continuously 
with prolonged incubation. In other words, carbon availability and 
lability did not certainly decease continuously. Such as, we can see 
that the parameter A for the sand fraction increased from the 
second to the third cycle both for grassland and forest soils 
(Figure 2). Tian et al. [71] also observed the phenomenon and 
affirmed that a pulse-dynamic pattern of soil respiration and 
increased carbon availability during a long incubation could exist. 
Overall, these results suggest that it is unreliable to judge whether 
the decomposition of recalcitrant and labile carbon have different 
temperature dependencies by comparing Qjo values after different 
incubation time. 

Conclusions 

This study presented three implications for understanding the 
soil-carbon response to climate warming. First, the Qjo of SOC 
decomposition in soil particles decreased with increasing particle 
size (clay>SLlt>sand) for grassland and forest soils. This result 
suggested that SOC detxmiposition from fme-textured soils rich in 
clay or silt content could be more sensitive to warming than coarse 
sandy soils. Second, Qjo decreased with increasing temperature, 
thereby implicating that the SOC might be more vulnerable in low 
MAT regions than in high MAT regions under fiiture warming. 
Third, the soil nitrogen availability facilitated a temperature 
dependence of SOC decomposition, suggesting that atmospheric 
nitrogen deposition could stimulate positive feedback of soil 
respiration on climate warming. 
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